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The renormaUzation of a scalar field theory with a quartic self-coupling (a A</>^ theory) via adi- 
abatic regularization in a general Robertson-Walker spacetime is discussed. The adiabatic coun- 
terterms are presented in a way that is most conducive to numerical computations. A variation 
of the adiabatic regularization method is presented which leads to analytic approximations for the 
energy-momentum tensor of the field and the quantum contribution to the effective mass of the 
mean field. Conservation of the energy-momentum tensor for the field is discussed and it is shown 
that the part of the energy-momentum tensor which depends only on the mean field is not con- 
served but the full renormalized energy-momentum tensor is conserved as expected and required by 
the semiclassical Einstein's equation. It is also shown that if the analytic approximations are used 
then the resulting approximate energy-momentum tensor is conserved. This allows a self-consistent 
backreaction calculation to be performed using the analytic approximations. The usefulness of the 
approximations is discussed. 



I. INTRODUCTION 



The study of free quantized fields in curved space has been a remarkably fruitful endeavor, par- 
ticularly in the applications that have been made to black hole and cosmological spacetimes [Q]. 
Much less has been done regarding interacting fields in curved space. However, interacting fields 
are very important since all real fields in nature appear to have interactions. Interactions also play 
an important role in cosmological models such as infiation, being required in some cases for the 
infiaton potential to have the right form and also for the thermalization that is necessary to reheat 
the universe after inflation. Interactions can also signiflcantly enhance the particle production that 
often occurs for free fields in curved space. 

The study of interacting quantum fields in Robertson- Walker (RW) spacetimes is of great impor- 
tance as well in understanding quantum fields in Minkowski space. It is well known that following 
a relativistic heavy ion collision, the quark-gluon plasma produced eventually undergoes a chiral 
phase transition. A good approximation to describing the dynamics of this system is provided by 
the linear a model, and by assuming that the expansion is mostly radial Written in terms of 

the spherical hydrodynamical fluid coordinates the system is equivalent to an interacting quantum 
scalar field (mean field plus fiuctuations) in a RW spacetime [Q. 

Perhaps the simplest interacting quantum field theory in four dimensions is a scalar field with 
a quartic self-coupling, often called the "A(/>*" theory. There is a long history of the study of this 
theory in curved space. The original investigations centered on renormalization. Drummond [^, 
Birrell and Ford j^Q, Bunch, Panangaden, and Parker and Bunch and Panangaden Q investi- 
gated the renormalization of the theory in various cosmological spacetimes using techniques such as 
dimensional regularization. Bunch and Parker showed that the theory is renormalizable in an 
arbitrary spacetime to second order in the coupling constant A. Birrell extended their work by 
using momentum space techniques and computing self-energy graphs to second order in A. 

Along with studies of the renormalization of the theory, various calculations have been undertaken. 
For example. Ford and Toms investigated phase transitions caused by one- loop radiative corrections 
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in an expanding universe jl^. The one-loop finite temperature eff'ective potential for a A(^'' theory 
in a RW universe was calculated by Hu, under the assumptions that the rate of change of quantum 
fluctuations is much greater than that of the mean field and the expansion rate of the universe |^ . 
Ringwald investigated the evolution of the expectation value of the quantum fiuctuation (i/)^) at 
one-loop order in a spatially flat RW universe [ [L4[ . 

The quantity (■i/)'^) can be used to determine the backreaction of the quantum fluctuations on 
the mean (or classical) scalar field as it appears as an effective mass for the mean field at one-loop 
order. However, to determine the backreaction of the scalar field on the spacetime geometry, one 
must compute the renormalized energy-momentum tensor for the field. The renormalization of the 
energy-momentum tensor to one-loop order in a spatially fiat RW spacetime was discussed by Paz 
and Mazzitelli They displayed the renormalization counterterms which were obtained using 

point splitting, adiabatic regularization, and dimensional regularization. The divergent counterterms 
were displayed in the format of dimensional regularization. Mazzitelli, Paz, and El Hasi used this 
formalism in a calculation relating to the evolution of the inflaton and the reheating after inflation 
in the new inflationary scenario [ p^ . 

There are two ways that so called "nonperturbative" effects are usually taken into account. One 
is the Hartree approximation which works for a single scalar fleld jl^. The other is the large 
N approximation where a single scalar fleld is replaced with A*' scalar flelds which are coupled 
via a quartic interaction, which is invariant under the group 0{N) and is thus often called the 
0{N) model. Mazzitelli and Paz considered the renormalization in both the Elartree and large A'^ 
approximations in an arbitrary background gravitational fleld, using point splitting techniques, and 
adiabatic and dimensional regularization [h8|. More recent work on the backreaction of a scalar 
fleld on the background geometry in a RW spacetime has been done for inflationary models by 
Boyanovsky, Cormier, de Vega, Holman, Kumar, Lee, Singh, and Srednicki |l9|-pl| and by Ramsey 
and Hu ||. 

In this paper we derive a set of renormalized equations that can be used to determine the evolution 
of a quartically coupled scalar fleld with arbitrary mass and curvature coupling to one-loop order in 
a RW spacetime. We also derive expressions for the unique components of the renormalized energy- 
momentum tensor that can be used to determine the backreaction of the fleld on the spacetime 
geometry. To renormalize we use the method of adiabatic regularization which is particularly useful 
for deriving a set of equations which are to be solved numerically. We display the counterterms 
for the energy-momentum tensor and the quantum contribution to the effective mass of the mean 
fleld (for arbitrary mass and curvature coupling). Previously the adiabatic counterterms for the 
energy-momentum tensor have been displayed by Paz and Mazzitelli but only in the context 
of dimensional regularization, which makes it difficult to use them for numerical computations, and 
by Ramsey and Hu for the minimally coupled case (if one takes the one- loop limit of their 1/N 
expansion) p2[ . 

We discuss in more detail than has previously been done the conservation of the energy-momentum 
tensor for the full system (mean field plus quantum fluctuations). We show that the natural division 
of this tensor into a "classical" and a "quantum" piece leads to neither piece being separately 
conserved. We also show explicitly that the full energy-momentum tensor is conserved. 

We present a variation on the method of adiabatic regularization which has been used by one of 
us to develop an analytic approximation for the energy-momentum tensor for a free scalar fleld in 
a RW spacetime |2^. We use this method to derive analytic approximations to both the energy- 
momentum tensor of the quantum fluctuation and to the effective mass of the mean fleld. If the 
analytic approximation is used in the equation for the mean fleld and if it is also used for the 
"quantum" energy-momentum tensor, then the resulting set of equations results in a conserved 
approximate energy-momentum tensor. Thus the analytic approximations can be used in lieu of 
the full renormalized expressions in the mean fleld and semiclassical backreaction equations. They 
are useful for the investigation of vacuum polarization effects, but not particle production since 
particle production is a nonlocal phenomenon. Nevertheless, the approximations give important 
information, particularly if one wishes to estimate the conditions under which the loop expansion 
breaks down. We discuss the validity of the approximation and argue that it is likely to be most 
useful for massless fields. 

It is well known that the problem of solving Einstein's equation in the presence of quantum matter 
is not an easy one. On the left-hand-side of Einstein's equation one needs higher order geometric 
tensors ('^'_ff^i/, '^'ff^^, and iJ,ji/) that involve up to fourth order time derivatives of the metric 
Qiiv 1^ • On the right-hand-side one encounters the expectation value of the energy-momentum tensor 
of the quantum field in a certain quantum state (for which there is no a priori rule to be determined 
or chosen), a quantity that is ultraviolet divergent. One has to regularize and renormalize (T^^) in 
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such a way that it remains covariantly conserved, as the left-hand-side of Einstein's equation is. This 
is why the problem of a fuU backreaction of the quantum matter on the spacetime geometry is so 
difficuh. One has to make use of regularization methods that are suited for a numeric computation 
(as the dynamical equation for the quantum field is most likely not to have analytic solutions), and 
at the same time guarantee the covariant conservation of the renormalized value of In this 

paper we use adiabatic subtraction to fulfill both requirements, and to set up all the formalism and 
techniques required to perform a backreaction calculation for an interacting theory (mean field plus 
fluctuations) in a general RW spacetime. 

Before involving ourselves in analysing the backreaction, it will be helpful to estimate how big 
or small the quantum efi^ects on the geometry and the mean field are. In this paper we introduce 
and describe the analytic approximation as such a tool. It yields a renormalized and covariantly 
conserved energy-momentum tensor for the quantum fluctuations, (that carries no information, 
whatsoever, about particle production effects), that can be evaluated in the spacetime geometry 
that is a solution of the Einstein's equation determined by the mean field. This paper presents all 
the technical details for such a calculation. Future work will consist of evaluating the analytically 
approximated energy-momentum tensor of the quantum fiuctuations in various scenarios, such as 
the reheating period of the infiationary regime of the early universe, and the spherical expansion of 
the quark-gluori plasma. 

In section |l| we introduce the conventions to be used, the background geometry, and derive 
the one-loop equations for the mean and quantum fields. We also compute the energy-momentum 
tensor of the system at one-loop, and show that it splits naturally in two terms: a "classical" 
and a "quantum" energy-momentum tensor. In section III we discuss the method of adiabatic 
regularization when the quantum fiuctuations have a time dependent mass, and derive the adiabatic 
order two and four counterterms that need to be subtracted from (V'^)ti and (TjuAu, respectively. We 
explicitly separate those new terms that were not present in the free case [p4|j2a] . In section |V| we 
introduce the analytic approximation as a way of estimating the importance of vacuum polarization 
effects, and as a first approximation to doing a full backreaction calculation. In section M we discuss 
the covariant conservation of the energy-momentum tensor. We show that the part of the energy- 
momentum tensor that depends only on the mean field is usually not conserved by itself, but that 
both the full energy-momentum tensor and its analytic approximation are covariantly conserved. 



II. BACKGROUND, CONVENTIONS, AND NOTATION 

We consider a quantum scalar field with self-interactions in a general RW spacetime. The metric 
of an RW spacetime can be written in the form^ 



ds^ = a^(ry) 



dr? - ~r dO, 

1 — ftr^ 



where rj is the conformal time coordinate, and «: = — 1, 0, +1 is the three-dimensional spatial curva- 
ture, corresponding to spatial Cauchy hypersurfaces that have negative, zero, and positive spatial 
curvature, respectively. 

The action of a scalar field with a quartic self-interaction is given by 



(2.1) 



(2.2) 



where g is the determinant of the metric, □ the D'Alembert wave operator given by □ = g^^V^V^, 
and R the scalar curvature of the RW spacetime. 

The equation of motion for the classical field (obtained by the principle of least action) is given 

by 

(^a + m^+^R+^^^^^^0. (2.3) 
The classical energy-momentum tensor is 



^Throughout this paper we use units such that h — c = 1. The metric signature is (H ) and the conventions for curvature 



tensors are R'^^s = ^'0-y,s — ■■■ and i?^^ = R'j^a,^. 
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If we quantize the theory $ becomes an operator. We then define the mean (or background) field 
by the equations 



where the expectation value is taken with respect to the initial state of the system (in the Heisenberg 
representation). 



Taking the expectation value of Eq. (2.3) and noting that 

(-l-^) = 4^' + 30^(7/') + 30(1/'') + (^^) = 0^ + 30(7/>') + {tp'') 
we find the following equation for the mean field 

(□ + m' + Ci?)0 + I (-^^ + 30(^') + (^')) = . 



In the same way, by subtracting equation (2.7) from equation (2.3), we obtain the equation of motion 
for the quantum fluctuation tp 

(□ + m' + ^R)tl; + ^ (3</)V " 30(1/'') + + tp'-" - (i/)^)) = . 

I f we trun cate at one-loop (free fleld theory for the quantum fluctuation ?/;), the equations of motion 
(E^ and (Ei) become 



3! 



{a + + ^R)iP + ^(j>^xP ^ . 

The expectation value of the energy-momentum tensor can be broken into a "classical" and a 
"quantum" part. The classical part is given by 



{T^,f = {l- 2£,)d^^d^ct> + (2^ - l/2)ff^,ac09' 

^2 A 



2e</'V^v„ 



while the quantum part is 

{T^,.f = (r„.)„ = (1 - 2^){d^M^iP) + (2C - 1/2)3^,. {9, V9"V) - 2C(^V„V,.V) 



Equations (2.9a- 2.10b) describe our system at one-loop order (mean field plus quantum fluctua- 
tions ■!/;, which contribute to the effective mass of 0). It is well known [p6lp4] that to make sense of 
these equations one needs to regularize the theory, that is, on the one hand, define a way to obtain 
from the bare parameters (m. A, and ^) the renormalized ones, and on the other hand, regularize 
the divergent quantities {tp^)u and {T^^)u, to obtain the physically finite energy-momentum tensor 
of the system. In the next section we discuss these issues. 



(2.4) 



(2.5a) 
(2.5b) 



(2.6) 
(2.7) 

(2.8) 

(2.9a) 
(2.9b) 



(2.10a) 



(2.10b) 



III. ADIABATIC REGULARIZATION 

Using dimensional regularization it has been shown that for a A0^ theory in a general spacetime, 
the bare and the renormalized parameters are related in the following way 

2 2 3Xr 2 

mg = rriR - —— —mji , 

87r-^(n — 4) 

6'^" 6 87r2(n-4) V*"* 6 

As = \r — „, ^ . . (3.1) 
87r^(n-4) 
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Here n is the number of dimensions the spacetime has been analytically continued to. Thus even for 
a general spacetime the counterterms for the renormalization of m^, ^, and A are constant in space 
and time as expected. 

In principle one would prefer to renormalize at the level of the effective action and then vary 
that action with respect to (j) and the metric to obtain the renormalized equations of motion 
and energy-momentum tensor, respectively . However the computation of the one- loop effective 
action for arbitrary mean fields (f> and arbitrary spacetime geometries is quite involved, and will not 
yield an intrinsically different answer from that obtained by looking at the one-loop field equations. 
For this reason, we consider a different method called adiabatic regularization p^-|3H which works at 
the level of the field equations. Another advantage of adiabatic regularization is that it is well suited 
to perform numerical calculations [ p2||3^ . In adiabatic regularization the divergences in quantities 
such as {ip'^)u and (^^v)u are computed using a WKB expansion for the modes of the quantized 
field V'- These terms are then subtracted from the unrenormalized (bare) expressions with the result 
that 

= (^'>" - {^'')a<i , (3.2a) 
{T^.)r = {T^u}u - {T^.)ad , (3.2b) 

where the subscripts u and ad stand for the unrenormalized (bare) and the adiabatic value, respec- 
tively of (?/>^) and (r^i^). This procedure has been shown to be equivalent to point splitti ng fo r free 
scalar fields in a RW spacetime |^,^. For the quartically coupled scalar field Eqs. (2.9a) and 
( 2.9b| ) then become 



(□ + ml + e«i?)0 + + = , (3.2c) 



(a + ml + ^RR)ij+^cl>^ij = . (3.2d) 

In what follows we only consider the renormalized values of the coupling constants m, ^, and A, so 
we drop the subscript R for these quantities. 

We assume that th e me an field is homogeneous (j) = (j}{ri), as a RW spacetime is homogeneous and 



isotropic. Then Eq. (2.9a) becomes 



,/ . 2 / 2 , ^ r) , ,2 , , , 2~ 



a 



+ 2-0' + a^(ml+ ^bR + + -^ir) ) = . (3.3a) 



Here primes denote deriva tives w ith respect to rj. The rj-q component of the classical renormalized 



energy-momentum tensor (2.10a) is given by 



= + + 3C f ^ + + + 4i^0^ , (3.3b) 



2^ ^ ^ a^^ ^ \ a'^ j ^ 2^ 4! 

and the trace is 

{T)% = (6C - 1)^<^V' + ^^H" + ^^-H' + ^R't>' + 2m^0^ + . (3.3c) 

a 3! 



To determine the renormalization counterterms used in adiabatic regularization, {{ip^)ad: {Tfj,„)ad), 
we first review canonical quantization in a RW spacetime. We then discuss the WBK expansion for 
the modes of the quantum fluctuation field ■;/), and compute the adiabatic counterterms needed to 
renormalize and {T^u)u- 

Since at the one-loop level the quantum field tp is a, free field with an effective mass of the form 
+ it can be expanded in the following manner Q 

V'(^) = ^ / dA(k) [flk Fk(x) Mv) + al (x) n (r,)] , (3.4) 

where the measure is given by 

/ d/i(k) = / d=*k if K = , 

= /o^°°dfc^,_^ if K = 

= y, , if K^+1 . 
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The spatial part of the mode function, YkCx), obeys the equation 

A(^'Fk(x) = -(fc2-K)Fk(x), 
and the time dependent part /fc(r?) is a solution to the equation 

+ ^a^<j>^ + - l/&)a^R = . 

Here uj^ = + m^a?, pri mes denote derivatives with respect to the conformal time r), and for 
spacetimes with the metric ( |2.l[ ) the scalar curvature is given by 



The unrenormalized (bare) expression for the quantum part of the energy-momentum tensor 
(|2.10b[) is HQ 



47r' 



/ dM(fe) j (m^a^ + \^<i?) \hf + 6 (C - i) 



/fcl (/fc/fc +/fc/fc) 



,2, 22, ^2,2," 1 , 

fc -l-ma + — atp H 5- + 

2 a 



(C-i)aMl/. 



Note that the pressure is not an independent quantity, and can be obtained from knowledge of the 
energy density and the trace. 

The equation for the mean field ( ^.9a| ) also contains the quantity {ip^). The unrenormalized 
expression for it can be written in terms of the mode functions as follows 



27r2a2 



In these expressions the measure is 



if K = 0, -1 
if K = +1 . 



To determine the adi abatic counterterms needed to renormalize these expectation values we solve 
the mode equation (3.6) using a WKB expansion. To obtain this expansion we first make the variable 
transformation 



exp 



dr]' Wkin) 



Substituting Eq. (^) into Eq. (|3J) yields 



2 1 / _ 3 
""21m/, 2 Wi 



This equation is then solved iteratively with uuk being of adiabatic order zero and the next two terms 
on the right hand side being of adiabatic order twon. Thus to second adiabatic order the solution 
for Wk is 



Wk=LOk + 



1 



T / ff 0/2 

1 I uJk 3 ujk 



2 \LUk 2 Lul 



(3.5) 



(3.6) 



(3.7) 



(3.8a) 



(3.8b) 



(3.8c) 



(3.9) 



(3.10) 



(3.11) 



^The \a^(f>^ term is considered to be of second adiabatic order because only terms with up to two time derivatives of (f> are 
needed to cancel divergences in {T^v)u |lq|. 
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To renormalize {tp'^)u we use a second order adiabatic expansion for the modes. A fourth order 
expansion is necessary to cancel the divergences in {Ty,i,)u [E4lj25l. The renormahzation counterterms 
for these quantities are 



dfiik) 



{T^ri)ad — {T^ri)ad + 



+ 



2 \ J,2„2„/2 



\{(l)-'a-' + 2(t)(t)'aa + K(t>^a'') Sm-^A^^a^a 



+ ■ 



(r)''''-(r)<>^+4^y g^33 + 32^(3A<^" + 4a + 4a ^ 

, o ,2 2 /2 I oj,2 3 // I TGAJ.' 3 "^^"^^/AJ.^ 4 /2 , ,2 5 // , nj, /' 5 /\ , 35m^'Xa^'a' 

+ Sep a a + Sep a a + Ibqxp a a ) — ^ ^ ^ (40 a a +000+ 2(?i0 a a ) + ■ 



/'OJ.J.' ' I J.2 " I J,/" 2 . //2 2 . 1,2 2\ 

--^^ — ^[Zq>q> aa + (p aa + (pq> a + (p a + nq) a ) 



(3<j!>'^a'^a''^ + 4(j)(j)' a^a + 2<jf>'^a'^a" + Kif'^a*) 



3m^A ^^^2^2^/2 , 3 / , o^2„3_// , J,2„4^ 15m'* Afli^a^'a' 



8^1 



The adiabatic counterterms for the free field are I p^jis) 



47r2 
1 



47r2a4 ^ 

128ajii 



d^(fc) "I <jjk + 



8a;5 32t^j^ 



3 / a 



-yZa a a ~a a + 4aa a 



12 It /4\ 



\iSiJl 



[aa a + a ) 



LOk \ a 



3m a' 3m f , ,,, ,,2 a' 



+ 



n ' I" n 

la a —a 



15m'' / /2 // , o /4 , 2 /2\ , 105m^a^a"' 
-=- [4aa a + 6a + a I + 



+ u 



9 f2a'a"' a"^ 4a' ^a" 2Ka 



+ + 



27m2 / a'^a" 



47r2a4 



dAt(fc) 



22 42 r64/2 42 

m a m a „ ,2-. 5m a a ma „„ ,./,„,„ /,2n 
H — : — + a ) = =r(aa + 4a a + ,ia ) 



4^i 



,,,1-1^, -;///, /2//,o2//2,o 

H Q-(4a a a + liiaa a + 6a a + Ja 

32Wfc 



8^1 



16^1 



,4. 231m* a*, ,2 „ ,4, 1155m^''a®a' 
Ti — {o.a a +a H rr. 



+ 



a / II I 

Ola a 



LOk \ a a 



3m2 



(2aa" — a'^ + Ka^) + 



n 4 2 /2 o 2 / rz 

am a a 6m I ,,,, a a 



+ 



+ 



15m /. 2 / /// I o 2 //2 I o 12 II /4 I 3 // I 2 /2\ 

=- i4a a a + 6a a + oaa a — a + na a + na a 

4loI ^ 



105m a 12 II , 14, , 2 /2x 
H — n — (Saa a + 5a + Ka a ) 



+ 



945m*'a'*a' 



9 /a"" 4a'a"' 3a"" 6a' "a" 2Ka" 2Ka' 



+ 



27m2 



/ 2 // 

, I III , 112 ^ a a / 2 22 

4a a + 6a — b h 4ftaa — Ina + a n 



135m / 12 II , 2 i2\ 
= — \ aa a + na a 



Equations ( |j.l2b| ), ( |3.12b[ ), and ( ^.12c| ) [together with ( ^.12d| ), ( |3.12e| ), and ( |3.12| )] give the adiabatic 
counterterms needed to obtain the renormalized expectation value of the quantum fluctuations and 
of the quantum energy-momentum tensor for an interacting quantum scalar field in a general RW 
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spacetime, for the case of a homogeneous mean field ijiiv)^- 



IV. ANALYTIC APPROXIMATIONS 



As it stands the method of adiabatic regularization can be used to compute the quantities {'iP^)r 
and {T^v)r- This then allows one to obtain a self-consistent solution to the mean field and mode 



equations (2.9a) and (2.9b) in a background RW spacetime (in the test field approximation) or 
to these equations plus the semiclassical backreaction equations in a RW spacetime. But before 
getting involved in the backreaction problem (and the difficulties this presents), it would be very 
useful to have a way of estimating how big or small the quantum corrections are beyond test field 
approximation. In this section we discuss this issue and present such an approximation. 

For a free quantum scalar field in a general RW spacetime, one of us has already shown that there is 
a way of defining a certain approximate energy-momentum tensor that is covariantly conserved |^ . 
In the present paper we extend this method to interacting quantum fields, and discuss its advantages 
and limitations. 

The analytic approximations result from a revision of the method of adiabatic regularization that 
simplifies the calculations in a RW spacetime with compact spatial sections, and in the process 
yields approximations for the quantities {iJP')r and {T^v)r. These analytic approximations give 
information about vacuum polarization effects but not particle production, since particle production 
is a nonlocal phenomenon. However they also make it possible to solve the mean field, mode, and 
semiclassical backreaction equations in an approximate manner, which goes beyond the test field 
approximation, and does not involve a full backreaction calculation. This can be useful, for example, 
if one wishes to determine under what conditions vacuum polarization effects will be important and 
what infiuence they may have on the mean field and the spacetime geometry. In particular, we believe 
that these analytic approximations may provide important information for reheating calculations, 
just before particle production from the infiaton field takes place. The analytic approximations 
provide a natural way to estimate the change in the vacuum polarization energy of the infiaton field. 
During the inflationary regime the mean fleld 4> dominates the energy density of the universe. It 
is reasonable to expect that when the vacuum polarization energy is of the order of the mean field 
energy density, the infiaton field will switch from the slow-roll regime to the oscillatory behavior, 
that will eventually lead to particle production. In this paper we restrict ourselves to presenting 
the analytic approximations and our method for obtaining them. In future work we will present the 
applications to reheating. 

In order to derive the analytic approximation, we first improve on the method of adiabatic regu- 
larization by expanding the renormalization counterterms in inverse powers of k, keeping only terms 
which are ultraviolet divergent. For the case of compact spatial sections (k = -1-1) the integral is also 



changed into a sum. We call the resulting expressions {'ip'^)d and (T^j^) 
RW spacetime they have the form 



respectively. In a general 



,2 / dM(fc) ^ 



+ 



dm - 



44 2\/,2 4 \2/,4 4 

m a m X(j> a X (j> a 



8 



8 



32 


2 








'2 \ 1 






)(^ 


-— 1 

a2 J 



(4.1a) 



(4.1b) 



^Mazzitelli and Paz presented these counterterms with the points separated and with the divergence structure given in the 
context of dimensional regularization |^ . Ramsey and Hu carried out the adiabatic expansion up to order four for ^ij = in 
the context of the leading expansion p^ . 
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2 2 / a 



A 



3m a 1" + T 



>H — 



II 2 , r±l'i 2 



+ 12(j)(j)' aa + ()<f?aa" + &ti(f)^ a 



(4.1c) 



with 



/ d^i(fe) = J^°°dkk'^ if K = 0, -1 



if «: = +! 



Here e is an arbitrary lower limit cutoff and 



(1) 



1 2 



(4.2a) 



In a RW spacetime it has the components 



(1) W 



a" 



+ ■ 



(1) 



a' 
J „iii 



+ 



+ 



a' 



+ 



+ 



+ ■ 



The renormalized e nergy -momentum tensor is then computed by subtracting and adding the quan- 
tity (r^^)d to Eq. (E2M, with the resuh that 



(4.2b) 
(4.2c) 



{Tj_iiy^ an = {Tf_iu') d — {Ti_ii/^ ad • 



(4.3a) 
(4.3b) 
(4.3c) 



In general {T^,^)„ must be computed numerically while {Tfj^)an can always be computed analytically. 
The result is 
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(4.4c) 



a" 1 
487r2a3 47r2 

a - l/6)i^ 
87r2 



m'' k{k + 1) 
"1^ 24^ ^ ~ 



log 



2e\ k{k + 1) 



(log e + 7) 



log 



f 2e\ ^(^+1) _^ ■ 
— n (loge + 7) 



V 



(4.5) 



9 



r] I an — 



28807r2 
647r2 

,2 



+ 



5 + 2log 



(f) 



+ K(K + l)(7 + loge) 



m^K{K — 1) 



1927r2a2 



2887r2 



^ «:(/^ - 1) 
327r2a4 



+ 



IBtt 



, G% ( 3 + 2 log f ^) + ^(^ + 1) (7 + log + ^"""^ 



327r2 V 



V2e 
2 



87r2a2 



21og(g)+.c(. + l)(7 + loge))-^ 



+ 



{T)L = 



28807r2 



+ 



6 

+ ^ 



m 



Sum 
9 

'87r2 



4a a' 



a" 



2887r2 
87r2a4 



2887r2 
1) / a 
a 



G 



967r2a2 



167r2 
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Here G^^ is the Einstein tensor with components 



G = 



3a' 



6a" 



3k 
a^ 



(4.6) 



(4.7) 

(4.8a) 
(4.8b) 



and Hu^ is the tensor 



with components 



3a' 



6Ka' 



-RpaR'"^ g^iv + T-R 5mi' 



3«;^ 



(3) 



H = 



12a' 



+ 



12Ka ' 



12Ka' 



(4.9a) 

(4.9b) 
(4.9c) 



For a massive field = m, while for a massless field is an arbitrary mass scale. However, in the 
massless case each of the terms containing log/i has as a coefficient a multiple of the tensor '^^'^H^,^, 
which comes from an term in the gravitational Lagrangian. Thus the terms containing log ^ 
simply correspond to a finite renormalization of the coefficient of the R^ term in the gravitational 
action. 

Note that if k = +1, {4''^)d and (Tp.v)d contain a sum over fc, while {'ip^)ad and {T^,y)ad contain 
an integral over k. Thus, either the integral must be converted to a sum or the sum to an integral. 
We have converted the sum to an integral using the Plana sum formula [p5M38|. This formula is 



^ fin) = i /(m) + 



n — m 



+ 00 /' + 00 

da; f{x) + i 



dt 



[f{m + it) ~ f{m - it)] 



Because of the way {Tp^)d is defined, the third term in the Plana sum formula can be computed 
exactly. In the traditional form of adiabatic regularization one would convert the integral in the 
adiabatic counterterms to a sum using the Plana sum formula and then substitute the result into 
Eq. (3.2b). However, if this is done then, for a massive field, it is not possible to compute the third 
term in the Plana sum formula analytically. Thus the computation of the renormalized energy- 
momentum tensor is simplified by our method in the n — +1 case. Clearly the same simplification 
would occur if one was using compact spatial sections for fi: = OorK = — 1 RW spacetimes. 



(4.10) 
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V. CONSERVATION OF THE ENERGY-MOMENTUM TENSOR 



In this section we show that the renormalized energy- momentum tensor for the A<^'* theory (at 
one-loop) in a RW spacetime is covariantly conserved. In a RW spacetime there is only one nontrivial 
conservation equation which is 

We have shown that the energy-momentum tensor for the field can be divided into a classical 
and a quant um p art. First consider the cl assi cal part whi ch in a RW spacetime has the components 
( 3.3b| ) and (3.3c). Substituting into Eq. (5.1) and using (3.3a) we find 



4a' „ a' A , , o, 

a a 2 



Thus the classical energy-momentum tensor is conserved only if there is no quantum correction to 
the effective mass of the mean field. 

The energy-momentum tensor for the quantum part consists of the difference between the un- 
renormalized part a nd the adiabatic counterterms. The components of the unrenormalized part are 
given in Eqs. (3.8a) and (3.8b). If they are substituted into Eq. (5.1) and the mode equation (3.6) 
is used then one finds 



If one substitute s the adiabatic counterterms (3.12b) and (3.12c) into Eq. (5.1), and compares the 
result with Eq. (3.2a) one finds 



4a' a' A 

{T'^ri)ad,-q + — {T'^ v) ad — — (T) ad = (/>(!>' {'4''^) ad 



Combining these results and using (3. 2a) and (3.2b), one finds that the total renormalized energy 



momentum tensor, classical plus quantum, is conserved. 

One further finds that if the analytic approximation is used in place of {'iP^)r in the equation 
for the mean field, and if it is used for the quantum energy- momentum tensor, then the analytic 
approximate energy-momentum tensor is conserved. This means that one can use the analytic 
approximation to define a consistent set of equations for the mean and the quantum fields, and 
to be the source in the right hand-side of Einstein's equation to solve a first approximation to the 
backreaction problem. 



Thu s (?/; )an and {T^^)an can be used in place of (t/; )_b and {T^,j)r in Eqs. ( 2.9a ), ( 2.10a ), and 
(2.10b) to obtain an analytic approximation for the system. For «: = 0, —1 RW spacetimes the terms 
being approximated contain the arbitrary constant e. This means the approximation is not unique 
unless the coefficients of the loge terms vanish. It is important to note that this is only true when 
using {'ip^)an and (r^y)an as an analytic approximation. The e dependent terms do not appear in 
the fully renormalized expressions of these quantities . 

From the DeWitt-Schwinger expansion [M it is known that for a free quantum field in the large 
mass limit {iP^)r and {Th,^)r have leading order terms proportional to l/ni^. Thus, the analytic 
approximations for these quantities are not good approximations in this limit. Previous numerical 
work indicates that the relevant condition is likely to be ma <C 1. The analytically approximated 
quantities are also local in the sense that they depend on the scale factor and its derivatives at a given 
time rj. Therefore, they cannot accurately describe particle production effects which are inherently 
nonlocal. However for massless fields they should allow one to estimate how important vacuum 
polarization effects are, and how they qualitatively effect the evolution of the system. 



VI. SUMMARY 



We have used adiabatic regularization to renormalize a scalar field theory with a quartic self- 
coupling of the form A<^* in an arbitrary RW spacetime. We have found that the energy-momentum 
tensor can be naturally split into two parts, a "classical" contribution (which corresponds to the 
energy-momentum tensor of a classical scalar field with a quartic interaction in a RW spacetime) 
and a "quantum" piece (which corresponds to the energy-momentum tensor of a free quantum scalar 
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field with tlio time dependent mass rn^ + ^^). We liave displayed the ronormalization counterterms 
for both the energy-momentum tensor and the contribution of the quantum fluctuations to the 
effective mass of the mean field at one-loop order. We have directly checked to see if the energy 
momentum tensor is covariantly conserved and found that while the entire tensor is conserved, its 
classical and quantum contributions are not separately conserved. 

By using a variant on the adiabatic rogularization method we have derived analytic approxima- 
tions for the energy- momentum tensor and the contribution of the quantum fluctuations to the 
effective mass of the mean field. We have shown that the approximate energy-momentum tensor 
is covariantly conserved. Thus the analytic approximations can be used in a self-consistent way to 
find approximate solutions to the mean field and backreaction equations. The approximations can 
provide a useful tool for learning about vacuum polarization effects for massless fields. However, 
they are not useful for massive fields in the large mass limit. They do not give any significant 
amount of information about particle production, which is a nonlocal phenomenon. 

The approximations could be of use in reheating calculations, in particular the transition from 
the slow-roll dynamics of the inflaton field to the oscillatory behavior around the minimum of the 
potential, (which is believed to produce particles of lighter masses) , and in the context of rclativistic 
heavy ion collision as a way of estimating the physical energy density and pressure of the vacuum and 
thermal excitations. The advantage of this approximation is that one obtains analytic expressions 
for the quantum piece of the energy-momentum tensor, without need to solve exactly the mode 
equation (which is the most difficult part to implement in numeric computations). In this way, 
it is relatively easy to study the backreaction problem of the full system (mean field, quantum 
fluctuations, and gravitational field). 

We believe that this program can be carried forward, and improved easily. We plan to extend the 
present approach and approximations to two- loop order, and to the expansion. We also plan to 
perform some specific calculations with applications to reheating and early density perturbations. 
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